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Green's functions in a holographically dual finite temperature field theory. The quasinor- 
mal modes are defined for gauge-invariant quantities which obey incoming-wave boundary 
condition at the horizon and Dirichlet condition at the boundary. As an application, we 
explicitly find poles of retarded correlation functions of i?-symmetry currents and the energy- 
momentum tensor in strongly coupled finite temperature M = 4 supersymmetric SU(N C ) 
Yang- Mills theory in the limit of large N c . 
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1. Introduction 

While the statistical origin of black hole entropy remains a subject of active research, one may 
wonder if the celebrated analogy [1] between the laws of black hole mechanics and the laws 
of thermodynamics can be generalized to non-equilibrium processes. Holographic AdS/CFT 
correspondence ( [2-4], see [5] for a review) provides a suitable arena for such a generalization. 
AdS/CFT conjecture asserts that string theories on certain asymptotically anti de Sitter 
spacetimes are dual to quantum field theories in lower dimension. Since the low-energy limit 
of string theory is described by the appropriate supergravity, problems in general relativity 
can be mapped to problems in the dual field theory. According to the duality, asymptotically 
AdS background spacetimes with event horizons are interpreted as thermal states in dual field 
theories. Correspondingly, small perturbations of a black hole or a black brane background 
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are interpreted as small deviations from thermodynamic equilibrium in a dual theory. This 
particular entry in the holographic dictionary can be made precise by considering quasinormal 
spectra of asymptotically AdS spacetimes. 

Quasinormal modes (see reviews [6] and references therein) are solutions to linearized 
equations obeyed by classical fluctuations of a gravitational background subject to specific 
boundary conditions. The choice of the boundary condition at the (future) horizon is dictated 
by the fact that classically horizons do not emit radiation. Thus out of two local solutions 
near the horizon typically representing waves incoming to the horizon and outgoing from it, 
one chooses the incoming waves only. This choice of the boundary condition has profound 
consequences, making the boundary value problem non-Hermitian, and the associated eigen- 
frequencies complex. This, however, is exactly what one expects in a holographically dual 
theory, where small deviations from thermal equilibrium are described by dispersion relations 
which correspond to non-zero damping [7]. Mathematically, these dispersion relations ap- 
pear as singularities of the retarded 1 Green's functions in the complex frequency plane. The 
connection between quasinormal spectrum of AdS black holes and singularities of thermal 
correlators in dual quantum field theories was first noted and explored for 2 + 1 dimensional 
BTZ black holes in [8]. It was pointed out later [9] that, even for higher-dimensional sys- 
tems, imposing Dirichlet boundary conditions for scalar perturbations at asymptotic infinity 
ensures that quasinormal frequencies coincide with the singularities of the retarded Green's 
function in a holographically dual theory. 

Quasinormal modes for electromagnetic and gravitational perturbations are physically 
more interesting than those for scalars because the corresponding fluctuations couple to con- 
served symmetry currents in the dual quantum field theory. However, the relation of these 
quasinormal modes to correlation functions of the dual theory is not immediate: for example, 
choosing Dirichlet boundary conditions for gauge-dependent quantities such as metric per- 
turbations would be rather unnatural. Thus, we address the following question in this paper: 
considering computation of a quasinormal spectrum as a purely general relativity problem 
(independent of the holographic duality) , what variables and boundary conditions should one 
use in order to ensure that the resulting spectrum coincides with the poles of the correlators 
in the dual quantum field theory? (For related discussions, see [10-12].) 

For vector and gravitational fluctuations, a convenient approach similar to the one used 
in cosmology [13] is to work with gauge-invariant combinations of the fluctuations. As an 
example, consider gravitational fluctuations h^ u of five-dimensional AdS-Schwarzschild back- 
ground with translationally invariant horizon. According to the gauge theory /gravity duality 
the fluctuation couples to the stress-energy tensor on the boundary [5], and thus we expect the 
quasinormal spectrum of h^ u to be related to the poles of the retarded two-point correlation 
function G^^p of the stress-energy tensor. As we discuss in Section the two-point function 
of the stress-energy tensor in the dual theory is a sum of three independent components 

Gn V>a p(k) = S^^Giiko^ 2 ) + Q^u, a /3G 2 (k ,k 2 ) + L^ Q/3 G 3 (ko,k 2 ) , 
1 Choosing outgoing waves at the horizon, one obtains advanced Green's functions in the dual theory. 
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where k is the four-momentum, and S^^, L^ U)a p are the appropriate projectors 

which provide three independent Lorentz index structures. The correlation function is there- 
fore completely determined by three scalar functions Gi, G2, G3. On the gravity side, one can 
form three gauge-invariant combinations of the components of in such a way that each of 
them corresponds to one of the functions Gi, G2, G3. Each of these gauge-invariant variables 
satisfies a second-order ordinary differential equation whose connection matrix essentially de- 
termines Gi, G2, G3. Dirichlet boundary condition imposed on each of the gauge- invariant 
variables ensures that their quasinormal spectrum corresponds to the poles of G\, G2, G3. 

The paper is organized as follows. In Section || we discuss general Lorentz index structure 
of thermal correlators of conserved currents and stress-energy tensor in relativistic quantum 
field theories. In Section |3| we propose a way to identify quasinormal frequencies of asymp- 
totically AdS spacetimes with poles of the corresponding retarded Green's functions in the 
holographically dual finite temperature field theory. In Section Q we provide a detailed ac- 
count of scalar, electromagnetic and gravitational quasinormal spectra for the five dimensional 
AdS-Schwarzschild background with translationally invariant horizon. Using the approach of 
gauge-invariant variables, we explicitly show that one can define quasinormal modes whose 
frequencies coincide with singularities of retarded Green's functions in the dual theory, when 
the latter are computed using the standard AdS/CFT prescription. In the low-energy limit 
we reproduce earlier results [14-16] on hydrodynamic properties of 3 + 1 dimensional N=A 
supersymmetric SU(N C ) Yang-Mills theory. In the more general case, we numerically com- 
pute the positions (in the complex frequency plane) of the singularities of retarded correlation 
functions of global i?-symmetry currents and energy-momentum tensor in strongly coupled 
AA=4 supersymmetric Yang-Mills theory in the limit of large N c . Some technical details 
appear in two appendices. 

2. Field theory correlators 

We start by discussing Lorentz index structure of retarded Green's functions of conserved 
currents and energy-momentum tensor in relativistic quantum field theories in infinite, flat 
D-dimensional Minkowski space. Translation and rotation invariance are assumed to be 
unbroken symmetries of the theory. We shall be interested in retarded Green's functions of 
conserved symmetry currents, 

C^(x -y) = -i 0(x°-y°)(Mx), My)}) (2.1) 

as well as of stress-energy tensor, 

G ma p{x-y) = -id(x -y )([T^(x),T a ^y)}) . (2.2) 

The expectation value is taken in a translation-invariant state, so that the expressions can be 
Fourier transformed: 

C^x -y) = f e ik ^C, u (k) (2.3) 
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and similarly for Gavapip ~ u)- Here k = (ko,k) is a D-dimensional momentum vector, 
kx = kox + kx, and the metric 7]^ u is taken to be mostly plus. Expectation values of all 
global conserved charges are assumed to vanish in the equilibrium state, in other words we 
consider systems without chemical potentials. Then CPT invariance of the equilibrium state 
implies that 

CW(fc) = C^(k) , (2.4) 
(k) . (2.5) 

In addition, correlation functions of stress-energy tensor satisfy 

G^apik) = G v ^ af j(k) = G^p a (k) (2.6) 

because of the symmetry of T^ v {x). Conservation of J^x) and T^ u (x) implies that the 
correlation functions may be defined so that they satisfy the following Ward identities 2 



kfC^ik) = , (2.7) 

k»G^ aP (k) = . (2.8) 

If in addition the theory possesses scale invariance, then the correlation functions of stress- 
energy tensors satisfy an extra Ward identity 

t v G^ aP {k) = . (2.9) 

Hermiticity of and combined with rotation invariance implies 

C ftv (k ,k) = C^(-k ,k)*, (2.10) 

G^apiko, k) = G^ af s(-k , k)* . (2.11) 

2.1 Conserved currents 



In vacuum, the Ward identity (2.7) implies that the current-current correlation function 
Cfj, v (k) is proportional to the projector onto conserved vectors, 

P„ = ^ - ^ , (2-12) 

where k 2 = —k 2 + k 2 . All components of C^ u (k) are thus determined by a single scalar 
function, 

C^(k) = P^U(k 2 ) . (2.13) 



2 One may choose to define the correlation functions in such a way that local (in position space) counter- 
terms appear on the right-hand side of the Ward identities. The correlation functions defined in this way will 
differ from C M „(fc) and G^,»5(fc) by analytic functions of ko and k. 
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If the expectation value is taken in a state that has only rotation symmetry (such as thermal 
equilibrium state in canonical ensemble), it is convenient to split the projector P^ u into 
transverse and longitudinal parts, 

P„u = P^ + P^, (2.14) 
where P^ v and P^ v are mutually orthogonal {Pa^ 11 P^p = 0) projectors defined as 

^o T o = 0, P£ = 0, P?j=8ij-^, (2.15) 
P^ = P^-PT. (2.16) 

They satisfy k^P^ v = k^P^ v = 0, and therefore any expression constructed out of P^ x , P^ x will 
automatically satisfy current-conservation constraint. Therefore in the rotation-invariant case 
the current-current correlation function is determined by two independent scalar functions, 

C^k) = P* U T (k , k 2 ) + I* U L (k , k 2 ) . (2.17) 

When n T = n L = n, this expression reduces to the Lorentz- invariant form fl2.13| ). Also, 
it is not difficult to see that due to rotation invariance the k —* limits (with ko fixed) of 
transverse and longitudinal self-energies coincide, 

lim n T (/c , k 2 ) = lim U L (k , k 2 ) . (2.18) 

fc^O fc^O 

As an example, consider a four-dimensional field theory at non-zero temperature. Without 
loss of generality one can take the spatial momentum oriented along the x 3 direction, so 
that kn = (—u,0,Q,q), with k 2 = —to 2 + q 2 . Then the components of the current-current 
correlation function are 



C x x xl (k) = C X 2 X 2 (k) = U T (co, q) , (2.19) 



as well as 



2 2 

Cu(k) = ^^U L (iO,q), C tx s(k) = ^-^n L (L0,q), C x s xi (k) = ^^U L (uJ,q) . (2.20) 
uj z —q z uj z —q z to z —q z 

For a system in stable thermodynamic equilibrium at temperature T, the low-energy (uj<^T, 
g<CT) behavior of H T and T1 L is universal and is described by effective hydrodynamic theory 
(see for example [17]). In the approximation of linearized hydrodynamics, H T (u>,q) is non- 
singular as a function of uj because it does not couple to charge density fluctuations. On the 
other hand, correlators which involve conserved charge density must exhibit a hydrodynamic 
singularity whose dispersion relation satisfies uj(q)^>0 as q^O. The longitudinal self-energy 
U L (to,q) has a simple pole at u>= — iD Q q 2 , where D Q is the diffusion constant of charge Q 
associated with current J^ix). 
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2.2 Stress-energy tensor 

In vacuum, the two-point correlation function of stress-energy tensor may be written as a 
sum of five terms allowed by the symmetries fl2.5|) , (|2.6| ), which are proportional to rj^rj^, 
(VmVvP + V/j.pVva), (il^ukak/3 + k^k u r] af3 ), (rj m k v kp + rjuak^kp + r]^pk u k a + r]„pk^k a ), and 
k^kykakp. Only two linear combinations of these terms are consistent with the Ward identity 
( |2.8[ ); they can be taken to be P iwPa.fi an d P^aPyfi + P^fiPva- A convenient way to write 
G^afi{k) is 

G ma p(k) = P^ u P a fi G B (k 2 ) + H^fi G s {k 2 ) , (2.21) 

where 

1 1 

Hfiu,afi = n [P/J-aPufi + PfifiPua) ~ JJ—^PfivPafi 

is a projector onto conserved traceless symmetric tensors, which is constructed to satisfy 
r]^ 1 'H ni/a/3 = 0. As a result, a scale-invariant theory must have Gs{k 2 ) = 0, and the correla- 
tion function takes a simple form 

G^, a fi(k) = H ma p G s (k 2 ) . (2.22) 

If the expectation value is taken in a state that has only rotation symmetry (such as thermal 
equilibrium state in the canonical ensemble), it is convenient to split H^ va p into mutually 
orthogonal projectors constructed out of PT V , P^v- One relevant combination is 

Sfiu,afi = r {P]iaPyfi + PftaPyfi + P^fiP^a + P^fiPua) ■ (2.23) 

It satisfies k^S^afi = 0, and also if v S wa ^ = because of the orthogonality of P T and P L . 
It is not difficult to find another independent combination with the same properties, 

Qnv,afi = JJ—^ {[P > - 2 )P^v P afi + JJ^P^v P Ifi ~ (P^u P afi + P^pPlfi^j ■ {2.24) 

The projectors S^^p and Q^ u ,afi square to themselves, and are orthogonal to each other 
(Sfiu, a fi V aX V Qxp,ar = 0). Therefore the projector H^^fi can be split as 

P^fj.u,afi ^fiu,afi Q[iu,afi ^^v^afi i 

where L^ a p = H^ a p - S^ a p - Q^afi is orthogonal to both S ma p and Q ma fi- Thus the 
correlation function of energy-momentum tensor in a scale-invariant theory can be written as 
a sum over three independent index structures, 

G^ a p{k) = S^ a pGi{k ,k 2 ) + Q fl u,afiG 2 {ko,k 2 ) + L^ a p G 3 (k ,k 2 ) . (2.25) 

It is not difficult to show that rotation invariance implies that k — > limits (with k$ fixed) 
of the three independent scalar functions must coincide, 

lim Gi(Jfeo, k 2 ) = lim G 2 {k , k 2 ) = lim G 3 {k 0l k 2 ) . (2.26) 

fc^O fc^O fc^O 
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In a scale non-invariant theory, two extra scalar functions are needed to specify Gavap{k). 
They multiply two independent linear combinations of P^P^p, P^P^pi an d (P^vPa/3 + 
Pjj^Pup); one possible choice is 



PluPlp + \{P T ^ + ) C T (k , k 2 ' 



+ *W/3 G i( fc o,fc 2 ) + Q^, a pG 2 (k ,k 2 ) + L^ a pG 3 (k ,k 2 ) . (2.27) 

When Ct=Cl=Gb and Gi=G 2 =G%=Gs, this expression reduces to the Lorentz-invariant 
form (OH). 



As an example, consider a four-dimensional field theory at non-zero temperature. Choos- 
ing momentum to be k^ = (— w, 0, 0, q) as above, one finds the following components of the 
correlation function. The correlations of transverse momentum density are determined by 
Gi(u,q), 



GW, tai (A;) = - 



G tx 1 ,x 1 x 3 (k) 



1 q z 
2^V 

1 ujq 



2u 2 -q 



Gi(iv,q), 
Gi(u,q), 



G x i x 3, x i x 3(k) = -— = — yGi(u, 



2uo 2 -q 2 



(2.28) 
(2.29) 
(2.30) 



The correlations of longitudinal momentum density, energy density, and diagonal stress are 
determined by G 2 (u,q), Ci,(u},q), Cx(oj,q). For example, 



G tt ,tt(k) 



1 



3{u 2 



G 



tt.tx 3 



(k) 



ujq" 



3 (u 2 -q 2 ) 2 



2G 2 (oj,q) + 3C L (oj,q) 
2G 2 {u,q)+3C L (uJ 



Gtt^x 1 



1 



2G 2 (cj,q) -3C L {u,q)-3C T (u,q) 



6 {q 2 — uj 2 ) 

The correlations of transverse stress are determined by Gs(uj, 

G x i x 2 :X i x 2(k) = ~G 3 (u,q). 



(2.31) 
(2.32) 
(2.33) 

(2.34) 



For a system in stable thermodynamic equilibrium at temperature T, the low-energy (oj<^T, 
q<^.T) behavior of G^vaB is universal and is described by effective hydrodynamic theory 
(see for example [17]). In the approximation of linearized hydrodynamics, G 3 (u, q) is non- 
singular as a function of uj because it does not couple to energy density or momentum density 
fluctuations. On the other hand, correlation functions which involve conserved densities ex- 
hibit hydrodynamic singularities whose dispersion relations satisfy co(q)^0 as q^O. Function 
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G\{uj,q) has a simple pole at uj= — i^ v Q 2 , where j v is the damping constant of the shear 
mode, proportional to shear viscosity. In a conformal theory (when Cl = Ct = 0), function 
G%(u, q) has simple poles at u> = ±v s q — iT s q 2 , where v s is the speed of sound, and T s is the 
damping constant of the sound mode, also proportional to shear viscosity. 



3. Quasinormal spectrum and holographic duality 
3.1 Thermal correlation functions 

We will be interested in small fluctuations of a black p-brane, 



d s 2 = (r) -f(r)dt 2 + ^{dx 1 ) 2 + b(r)dr 2 . (3.1) 



i=l 



Metrics of this form arise as a result of dimensional reduction of higher dimensional super- 
gravity backgrounds. In addition, these backgrounds have non-zero values of various "matter" 
fields which we collectively denote by <f>^ suppressing all indices. Holographically dual theory 



is defined on the boundary (r^co) of (3.1), which is a flat p+1 dimensional Minkowski space. 

Translation invariance on the boundary implies that all fluctuating fields can be taken 
to be proportional to e - lujt + l i x ^ thus linearized fluctuations 5g^ u , 5<p of the background will 
obey a system of second-order linear ordinary differential equations. Generically, the system 
will be redundant, reflecting the gauge freedom (such as linearized diffeomorphisms) enjoyed 
by the fluctuation fields. Instead of fixing a particular gauge, we consider gauge-invariant 
combinations of the fluctuation fields. Let be these gauge-invariant variables which are 
constructed as linear combinations of fluctuating fields and their derivatives, excluding r- 
derivatives. Variables Zf. will obey a system of coupled second order linear ordinary differential 
equations (ODEs) which can in principle be diagonalized. Let Z(r) be such a gauge-invariant 
variable satisfying a second-order ODE. A local solution of the ODE near the horizon will 
generally be a superposition of incoming and outgoing waves. Since classically the horizon 
does not radiate, we choose the incoming wave boundary condition there. 

The solution obeying the incoming wave boundary condition at the horizon can be written 
in the basis of two local solutions at the boundary 3 as 

Z(r)=A<pi(r)+Bcp 2 (r), (3.2) 

where A, B are the connection coefficients of the corresponding ODE. Coefficients .A, B 
typically depend on the parameters (such as frequency and momentum) which enter the 
differential equation for Z(r). Near the boundary, the solution ( |3.2| ) becomes 

Z(r)=Ar- A - (l + ---)+Br- A + (! + ■■•), (3-3) 



3 By local solutions we mean the solutions obtained as power series around the corresponding singular point 
(r—oo in this case) of the differential equation. See for example [18] for a general discussion. 
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where A + , A_ are exponents of the ODE at r = oo, and ellipses denote higher powers of 
r. We consider the situation when the exponents are not equal, A + > A_, and A + can be 
taken positive. 4 

The action of the system expanded to quadratic order in fluctuations <5flw, ^ can be 
rewritten in terms of the gauge-invariant variables. On shell, the part of the action quadratic 
in fluctuations will reduce to the boundary term 



~ lim / dwdPq F(u, q) Z'(r) Z(r) + contact terms , 

r— >oo J 



(3.4) 



where F(uj, q) depends on the details of the action, and "contact terms" do not contain Z'(r). 

In holographic AdS/CFT duality, fluctuation 5(f) couples to a particular operator O of the 
dual theory at the boundary. Applying the Lorentzian AdS/CFT prescription [9,20] to the 
action (|3.4|) to compute the retarded correlator, and remembering that Z(r) is a functional 
of 5(f), we find 

B 

(OO)r ~ — + contact terms . (3-5) 

The poles of the retarded correlator correspond to zeros of the connection coefficient A. On 
the other hand, setting A = in Eq. ( |3.2| ) corresponds to a particular choice of boundary 
conditions for the fluctuation Z(r). From the general relativity point of view, this choice 
determines the quasinormal spectrum of Z(r). In other words, equation A = defines 
quasinormal spectrum for gauge-invariant perturbations which has the interpretation of the 
poles of retarded correlators in a holographically dual theory. 5 This argument is similar to 
the one given in [9, 22] in the case of a scalar fluctuation, the difference being the use of 
gauge-invariant variables in the present discussion. 



4 In AdS/CFT correspondence, A+ is equal to the conformal dimension of the operator that couples to bulk 
fields contained in Z. See also [19] for exceptional cases. 



5 One should add at once that in addition to zeros of A, singularities of the correlator (3.5) may also come 
from singularities of B. However, singularities of B are completely determined by the singularities of the local 
Frobenius solution <pi considered as a function of parameter(s) of the differential equation. Indeed, a general 
theorem [21] guarantees smoothness of a solution of a differential equation with respect to a parameter, if the 
equation and the boundary conditions depend smoothly on the parameter. Thus the solution Z is smooth, and 
singularities of B are destined to cancel the singularities of the coefficients of the series expansion in the local 
Frobenius solution tpi with respect to a parameter. Therefore, singularities of B are essentially determined by 
the recursion relations for the coefficients of the series expansion which defines tp\ . Let us illustrate this point 
using hypergeometric equation as an example. On the interval z £ [0, 1], the solution 2F1 (a, b; c; z) defined by 
its recursion relations at z — is related to two local solutions at z = 1 by 

2F1 (a, b;c;z) = A 2P1 (a, b; a + b — c + 1; 1 — z) 

+ B (1 - z) c - a ~ b 2F1 (c - a, c - 6; c - a - b + 1; 1 - z) , (3.6) 

where the connection matrix coefficients are given by 

r(c)r(c-a-6) r(c)r(a + b-c) 

r( c - o)r( c - b) ' r(a)r(6) 



To be in agreement with the scenario set by Eq. (3.2), let us further assume that c — a — b > 0. The poles of 



the correlator (p.5|) come from the poles of T(c — a) and F(c — b) (also corresponding to the Dirichlet condition 
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3.2 Black brane fluctuations 

According to the dictionary of the AdS/CFT correspondence [5], global symmetry currents of 
the dual field theory have as their sources boundary values of the gauge field on the higher- 
dimensional background (|3.1| ). Similarly, stress-energy tensor of the dual theory is sourced 
by gravitational fluctuations hav of the black brane. We shall be interested in quasinormal 
spectra of these fluctuations. We take the fluctuations to be of the form A^(r)e~ lujt+iqz , 
h^ u {r)e~' LUlt+iqz where z = x p . The fluctuations can be classified according to their transfor- 
mation properties under the "remaining" world- volume symmetry group 0(p— 1) acting on 

1 c— 1 

CO cc 

Let us restrict ourselves to the simplest case when the considered fluctuations do not 
couple to fluctuations of other background fields. Components At, A z do not transform 
under 0(p— 1), while components A a , a = x , .., x p ~ l transform as vectors. One can therefore 
distinguish two symmetry channels for electromagnetic fluctuations 

Spin (diffusive channel): A t , A z , A r (3.7a) 
Spin 1 (transverse channel): A a . (3.7b) 

A similar classification can be adopted for metric fluctuations. Components hu, ht z , h zz , h rr , 
htr, and h zr do not transform under 0(p— 1), components ht a , h za , and h ra transform as 
vectors, while h a a transform as rank-2 tensors. The tensor representation is reducible, for 
a symmetric h a p can be decomposed into the trace part 5 a ph/{p— 1), where h = J2 a ^aa (a 
singlet) and the symmetric traceless part h a » — 6 a ah/{p—l). Thus we have three symmetry 
channels for gravitational fluctuations: 6 

Spin (sound channel): htt, ht z , h zz , h, h rr , ht r , h zr (3.8a) 
Spin 1 (shear channel): hf a , h za , h ra (3.8b) 
Spin 2 (scalar channel): h a p — 5 a ph/{p— 1) . (3.8c) 

The 0(p—l) symmetry guarantees that equations for fluctuations belonging to different sym- 
metry channels decouple. Classification presented here mirrors the classification of the cor- 
relators in Section ||: diffusive and transverse channels of the U(l) fluctuation correspond 

A — 0), and from the poles of V(a + b — c). The latter are determined by the local solution near 2 = 1, 

2Fi(a,b;a + b- c+l;l- z) = 1 + - ^— (1 - z) 

1 + a + b — c 

afe(a +!)(&+!) _ a 

2(l + a + b-c)(2 + a + b~c) { ' 

the coefficients of which have poles at a + b — c = — n, where n is a positive integer. These are precisely the 



poles of the correlator (3.5) coming from £>. By setting a + b — c = —n + e and taking the limit e — > 0, one can 



show that the right hand side of Eq. (3.6) is in fact a smooth function of the parameters a, b, c. 

6 The name "scalar" for spin-2 fluctuations reflects the fact that the corresponding wave equation coincides 
with that of the minimally coupled massless scalar [23] . The names "shear" and "sound" (as well as "diffusive" 
above) reflect physical interpretation of the lowest quasinormal frequency (for a given symmetry channel) in 
the dual field theory. This will be seen explicitly in the next section. 
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to functions Tl L and II T of the current-current correlator. When the dual theory possesses 
conformal invariance, shear, sound and scalar channels of the gravitational fluctuation are 
related to functions G\, G2, G3 of the stress-energy tensor correlator. 

3.3 Gauge-invariant variables 

We now define gauge- invariant variables corresponding to classes (p.7|), (|3.8[). Combinations 



of the gauge field components invariant under the transformation — > — d^X include 
components of the electric field 

Diffusive channel: E z = qAt + ivA z , (3-9) 
Transverse channel: E a = LoA a . (3.10) 

Perturbations transform under infinitesimal diffeomorphisms as h^—^h^ — £, V — V V 
where = !;^(r)e~ lu ' t+igz are gauge functions, and covariant derivatives are taken with re- 
spect to the background metric ( |3.1|) . One may define the following gauge-invariant combi- 
nations linear in perturbations: 

Shear channel: Z\ = qH tx i + loH zx i (3-H) 

Sound channel: Z 2 = q 2 fH tt + 2uqH u + cu 2 H zz + q 2 f(l + ^-- H (3.12) 

V a f if J 

Scalar channel: Z3 = H x i x 2 , (3.13) 

where H tt = h a /af, H tz = h tz /a, Hij = hij/a / t), H = h/(p-l)a. From Einstein 
equations obeyed by the fluctuations, one obtains three independent second-order ODEs 
satisfied by Z%, Z2, and Z3. 

Quasinormal modes are defined as solutions to the second-order differential equations 
satisfied by the gauge- invariant variables E z , E a , Z%, Z2, Z3 obeying incoming wave boundary 
condition at the horizon and Dirichlet condition at the boundary. According to the above 
discussion, the spectra of complex eigenfrequencies obtained in solving the boundary value 
problem for fluctuations E z , E a , Z\, Z2, Z% coincide with poles of the functions T1 L , T1 T , G\, 
G2, G3, respectively. 

4. Thermal correlators in strongly coupled J\f = 4 SYM 

We now apply the approach outlined in Section || to the near-horizon limit of the non-extremal 
gravitational background of type IIB low energy string theory describing iV c parallel black 
three-branes. The background is given by the metric 

r 2 i? 2 
ds 2 = (-f(r)dt 2 + dx 2 + dy 2 + dz 2 ) + -^j^dr 2 + R 2 dQ 2 5 , (4.1) 

4/4 



where R is a constant, which depends on the number of D3 branes, R oc N c , and f(r 

r 4 




1 — rp/r 4 . The parameter of non-extremality rg specifies the location of the horizon, whose 
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Hawking temperature is T = r^jixR 2 . Upon dimensional reduction on S 5 , gravitational 



perturbations will effectively propagate on the five-dimensional part of the background (4.1). 
Introducing new coordinate u = r^/r 2 , the metric can be written as 

ds 2 = (-f(u)dt 2 + dx 2 + dy 2 + dz 2 ) + - — TTTT^rdu 2 , (4.2) 

u iu z j(u) 

where f(u) = 1 — u 2 . In these coordinates, the horizon is located at u=l, and the boundary 
is at u=0. In addition, the background is specified by the value of the self-dual five form field 

2(irTR) 4 



-(l + *)dt Adx Ady Adz A du , (4.3) 



with all other fields vanishing. The dual quantum field theory is N = 4 SU(N C ) supersym- 
metric Yang-Mills theory in 3 + 1 dimensional Minkowski space at large N c and large 't Hooft 
coupling [2]. The field theory is taken in a thermal equilibrium state at a temperature equal 
to the Hawking temperature of the background. Real-time thermal correlators of the con- 
served i?-symmetry currents and stress-energy tensor in this theory were considered in the 
AdS/CFT approach in [11, 15, 16]. Here we show that reformulating the problem in terms 
of gauge-invariant variables allows one to compute functions H L , H T , G\, G2, G3 directly 
by solving the second-order ODEs associated with each of them. In particular, quasinor- 
mal spectra of the gauge-invariant fluctuations determine poles of these functions in complex 
frequency plane. 

4.1 i?-current correlators 

Correlators of i?-currents in strongly coupled M = 4 SYM at zero temperature were computed 
in [24,25] by using the AdS/CFT correspondence. Similar approach can be taken at non-zero 
temperature [15]. On the gravity side of AdS/CFT, one considers an effective U(l) field 8 in 
the five-dimensional asymptotically AdS part of the background (4.2). This five-dimensional 



Maxwell field is essentially a graviphoton of the dimensional reduction of ( f4.1| ) on 5 5 . 

According to the discussion in Section thermal current-current correlators are deter- 
mined by two independent scalar functions, H T (io,q) and Tl L (u},q). Correspondingly, we 
expect the dual five-dimensional Maxwell system to reduce to two independent equations for 
gauge-invariant variables whose quasinormal frequencies determine the poles of the correla- 
tors. One of the quasinormal frequencies should be purely imaginary (at least in the regime 
q/T <C 1), reflecting diffusive relaxation of large-scale charge density fluctuations around 
thermal equilibrium state in the dual field theory. 

Maxwell's equations for the U(l) field are simply d A (^—g g AC g BD F CD )=0, where F CD = 
d c A D — d D A c , capital Latin indices run over t,x,y,z,u, and components of the five dimen- 
sional metric g AB are given by (|4.2[). Translation invariance for the t, x, y, z coordinates implies 



r The problem of additional singularities discussed in footnote || does not arise here. 

8 The _R-charges of the theory Q a = f d 3 xjo(x) generate global 5'J7(4) symmetry group, with a = 1 ... 15. 
In an equilibrium state without chemical potentials for the i?-charges, the correlation function of R currents 
j 1 ^ has the form C^ v {x) = 5 ab C IJ , v (x). The expressions of this section refer to C^ v {x). 
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that vector potential can be Fourier transformed, 

7i 



A c (u, t,x) = J pe iW+fa 4 c ( U) k) . (4.4) 

Choosing k = (— u, 0, 0, q), one can derive the following equations satisfied by transverse and 
longitudinal electric fields: 



K + -J E 'a + 1 n^ L E a = 0, a = x,y, (4.5a) 



E " + f( f - m E * + ^ *t S E * = < ( 4 - 5b ) 

f(vo 2 - q 2 f) uf 2 
where E a = voA a , E z = qA t + voA z , dimensionless parameters t» and q are defined as 

and prime denotes the derivative with respect to u. 

Gauge-gravity duality implies that all information about two-point i?-current correlation 
functions in the dual jV=4 SYM theory (in the large N c and large 't Hooft coupling limit) is 
contained in the solutions to the differential equations (4.5). For both equations, the singu- 
larity at u = 1 (the horizon) has exponents ±zt»/2 corresponding to the outgoing/incoming 
waves. To compute the retarded correlators, one has to impose the incoming wave boundary 
condition at the horizon [9] thus choosing —ixo/2 as the correct exponent. At the boundary 
(u = 0) the exponents for both equations are and 1, and thus solutions to equations (4.5) 
which satisfy incoming-wave boundary condition at the horizon behave near u = as 

E a (u) = -4( a ) (to, q) H h£( a )(»,q)uH , (4.7a) 

E z (u) = A {z) {to,q) + --- + B {z) {w,c,)u+-- - . (4.7b) 

The boundary action of the Maxwell system in the gauge A u = is 9 

S = lim ^ J jgf [A' t (u, k)A t (u, -k) - /(«) A'(u, k)A(u, -k)] (4.8) 
Using Maxwell's equations, the action can be written in terms of gauge-invariant variables as 



S = lim 



N 2 T 2 f duodq f f(u) 



u^o 16 J (2vr) 2 Wf(u) - to 2 
/(«) 



E z (u, k)E z (u, —k) 



., yE' x (u, k)E x (u, —k) + E' y (u, k)E y (u, —k)) + contact terms , (4.9) 



where "contact terms" do not contain derivatives of the electric fields. In order to find the 
correlation functions, one has to express the derivatives of the fields in terms of the boundary 



formalization of the five-dimensional action S = l/4g% J y^g F AB F AB is fixed by g% = 16tt 2 R/N^ [24] 
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values of the fields A^(k)=A^(u^0, k), E®(k)=En(u— >0, k) using the solutions (4/7); applying 
Lorentzian AdS / CFT prescription [9] , one finds 10 

5 2 S w 2 5 2 S N 2 T 2 B 



C aa {u),q) 



(a) 



6A&{k)6A%(-k) 5E£{k)5E£(-k) 



8A 



(a) 



Similarly, 



n 2 tWb (z) 



C zz (u,q) 



8(q 2 



m 2 )A (z) ' 

N 2 T 2 wqB {z) 



N 2 T 2 w 2 B 



8(q 2 

(*) 



vo 2 )A iz) 



8(q 2 



to 



OA*) 



Comparing these expressions with the general result ( |2,20| ), one finds 

N 2 T 2 B {a) (iO,q) N 2 T 2 B (z) (iO,q) 



V- T (LO,q) 



U L (u,q) 



(4.10) 

(4.11a) 
(4.11b) 
(4.11c) 

(4.12) 



8A( a ){u,q) ' ' ' ' 8A( z ){u,q) 

Thus the correlation functions are completely determined by the ratios of the connection 
coefficients of differential equations (4.5). In particular, poles of the correlators correspond 
to zeros of the coefficients A( a )(vo, q) and Ar z )(to, q). To find the zeros, we impose Dirichlet 
boundary conditions on electric fields at u = for the solutions to equations (4J3) which 
satisfy the incoming wave conditions at the horizon. Physically, the horizon acts as a perfectly 
absorbing surface, while the boundary acts as a perfect conductor. 

In order to determine the self-energies II T (a;, q), n^cj, q), one needs to know the solution 
to equations ( ]4 . 5| ) . Analytic solution is unknown, except for a special case of q = 0, when the 
Dirichlet boundary value problem can be reformulated as a problem of solving a transcen- 
dental algebraic continued fraction equation [11,26]. For q = 0, gauge- invariant variables Ei, 
i = x,y, z obey the same equation 11 



f uf 2 



0. 



(4.13) 



The solution to the Dirichlet boundary value problem for Eq. ( [4.13 ) can be found exactly, 
and is given by Heun polynomials [11]. The quasinormal spectrum is 



0, 



n(l 



n 



0,1,2,.... 



(4.14) 



For q / 0, quasinormal spectra of perturbations E a and E z (and, correspondingly, the poles of 
H T (ui,q) and II (oj, q)) differ from each other, and can be found numerically, as explained in 
Appendix [B|. A typical arrangement of quasinormal frequencies is shown in Fig. [I] (the poles of 
H L (uj,q) were found numerically in [11]). Quasinormal frequencies are located symmetrically 
with respect to the imaginary axis, as is expected from the singularities of the corresponding 
correlation function in the dual field theory, see equation (|2.10| ). 
10 See [9] for the definition of functional derivative in this context. 

n As they should, in accord with rotation invariance in the dual field theory, see Eq. (2.18). 
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Figure 1: Quasinormal spectrum of electric field fluctuations in the plane of complex to = oj/2ttT, 
shown for spatial momentum q = qjlixT = 1. Quasinormal frequencies on the left are defined by 



equation ( 4.5a ) for the transverse electric field, and coincide with the poles of H T (uj,q), as explained 
in the text. Quasinormal frequencies on the right are defined by equation ( 4.5b ) for the longitudinal 
electric field, and coincide with the poles of U l (uj, q). As q decreases, all poles stay at a finite distance 
from the real axis, except for the one marked with a large dot. This pole is purely imaginary and 
approaches the origin in the limit q — + 0. The presence of this special quasinormal frequency is 
a manifestation of the diffusive relaxation of i?-charge density fluctuations in the dual J\f=4 SYM 
theory. 



Hydrodynamic approximation 

In the hydrodynamic limit one can find analytic solutions to equations ( |4.5| ) as a series in 
to <C 1, q <C 1. Assuming first that to, q are of the same order, we get 



E a (u) =C a f{u)~ lX °l 2 
E z (u)=C z f(u)-^ 2 



1 + ito log 
1 + ito log 



1+u 
2 

1 + u 



+ ^: ( i_ n) + 0(n ,2 )C ^ 



(4.15a) 
(4.15b) 



where C a , C z are normalization constants. Imposing Dirichlet boundary conditions at u = 0, 
one finds that the equation E a (0) = has no solution compatible with the assumption to -C 1. 
Holographic interpretation of this fact is that the function U T (uj,q) has no singularities in 
the hydrodynamic regime. The condition E z (0) = leads to to = — iq 2 + 0(q 3 ). This is the 
lowest hydrodynamic quasinormal frequency of the Dirichlet boundary value problem for E z . 
Accordingly, in the hydrodynamic regime function n^u;,^) has a pole at 12 



-iD n q 2 



(4.16) 



where D R = 1/2ttT. Physically, it corresponds to the R-charge diffusion with the diffusion 
constant D R . By comparing the solutions ( 4.15| ) to the definition of connection coefficients 

12 One may question the validity of this result since it implies that t» ~ q 2 , whereas the solution (4.15b) 
was obtained under the assumption that to and q are of the same order. To check our result, we introduce a 
new parameter /i = t»/q ~ q. Solving Eq. (4.5b) perturbatively in < 1, (| C 1, where /i and q are of the 
same order, we find E z (u) = C z /-*^/ 2 [1 + iq (1 - uj/fj, + O(p)}. The condition E z (0) = a gain gives the 



dispersion relation vo 
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(4.7), one finds Ai a \ = 1, £>( Q ) = ito, Au\ = 1 + iq 2 /t», = ito — iq 2 /to. Substituting these 



connection coefficients into AdS/CFT results ( 4.10| ) and ( |4.11| ), one reproduces the i?-current 



correlators in the low-frequency approximation found earlier in [15]. 
4.2 Stress-energy tensor correlators 

Correlation functions of stress-energy tensor in the finite temperature N = 4 SYM at large 
N c and strong coupling were studied in [11,15,16]. Here we use the gauge-invariant variables 
approach to identify the correct associated boundary value problem, and to obtain new results 
for the correlators in the sound wave channel. 

To compute correlators of the stress-energy tensor in AdS/CFT correspondence, one 
considers metric fluctuations h AB of the supergravity background. To linear order in h AB , the 
Einstein equations are 

n ™ = ~W Hab ' (417) 

where TZa B is the linearized Ricci tensor evaluated in the background (14.21) . Translation invari- 



ance for the t, x, y, z coordinates implies that metric perturbations can be Fourier transformed, 
and classified according to their transformations with respect to the rotation group 0(2), as 
discussed in Section |3.2| . As we shall see shortly, quasinormal spectra of the three gauge- 
invariant variables ( |3.11| ) - ( [3. 13 ) appear correspondingly as the poles of the three functions 



Gi, G2, G3 in Eq. (|2.25[) which determines the two-point correlation function of stress-energy 
tensor in a scale-invariant theory. 

4.2.1 Scalar channel 

According to the discussion in Section [| the equation satisfied by the component h xy of 
the perturbed metric decouples from the rest of Einstein equations. The gauge-invariant 
function Z 3 = H xy = hy satisfies the equation for a minimally coupled massless scalar in the 
background (fl^), 

1+u 2 ro 2_ q2/ 

The exponents of this equation near u=0 are and 2, therefore asymptotic behavior of Z3 
near the boundary is 

Z 3 (u) = A {3) (1 + • • • ) + B (3) u 2 + • • • , (4.19) 

where ellipses denote higher powers of u. The relevant part of the boundary gravitational 
action (Eq.(6.9) of [15]) can be written using our notations as 

vr 2 7V r 2 T 4 , f dudq f(u) „, , T « , , 
S = limj ^l^z'^k^-k). (4.20) 



Z »-i^Lz; + ^l^ Z3 =0. (4.18) 



Proceeding as in Section 4.1 , we obtain the correlator 13 



G X y,xy — . • (4-21) 

2^.(3) 



3 Terms analytic in t» and q are ignored, even if they are divergent as u^O. 



- 16 - 



-4 


-2 




2 


4 










Ret» 






-1 


• 








-2 








• 












-3 










-4 






• 














-5 






■ 






ImtO 





Figure 2: Quasinormal spectrum of gravitational fluctuations in the scalar channel, shown in the 
plane of complex to = uj/2ttT, for spatial momentum q = q/2nT = 1. The quasinormal frequencies 
coincide with poles of Gs(uj,q), as explained in the text. As q — > 0, all poles stay at a finite distance 
from the real axis, as one expects from the absence of hydrodynamic singularities in Gs(u>, q). 



Comparing this expression to the general result ( 2.34 ), we nn d 



G 3 {uJ,q) = — r . (4.22) 

-4(3) 

Turning now to the connection between quasinormal spectrum and AdS/CFT correlators, we 
see that the condition ^4( 3 )(i»,q) = determines the singularities of G3. From the definition 
of connection coefficients Q4.19 ) it is evident that the condition „4( 3 )(t»,q) = is equivalent 



to imposing the Dirichlet boundary condition on fluctuations, ^3(^=0) = 0. 

The quasinormal spectrum of fluctuations obeying incoming wave boundary condition at 
the horizon and Dirichlet condition Zz(u=Q) = at the boundary was numerically computed 
in [11,26]. The spectrum o> n is discrete, (presumably) infinite, and almost equidistant. For 
q = its asymptotics for higher modes is well approximated by a simple formula 14 [26] 

oj n = 2vrTn (±1 - *) , n -» 00 . (4.23) 

A typical arrangement of quasinormal frequencies is shown in Fig. ||. Quasinormal frequen- 
cies are located symmetrically with respect to the imaginary axis, as is expected from the 
singularities of the corresponding correlation function in the dual field theory, see Eq. ( p. 11 ). 



Hydrodynamic approximation 

In the hydrodynamic limit xo <C 1, q <C 1, analytic solution to Eq. ( |4,18| ) can be found, 

Z 3 (u) = C7 3 f{u)-™/ 2 [1 + O(vo\ q 2 )] , (4.24) 

where C3 is a normalization constant. By comparing this solution with the definition of 
connection coefficients ( |4 . 1 9| ) , one finds A(%\ = 1 + 0(xx> 2 ), = ixo/2 + 0(xx> 2 , q 2 , toq). 



For an analytic approach to the asymptotic behavior (J4.23|) see [27] 
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Substituting these connection coefficients into the AdS/CFT result ( 4.21Q , one reproduces 
the correlator of transverse components of stress tensor in the low-frequency approximation 
found earlier in [15]. The equation »4(3)(t»,q) = does not have solutions compatible with 
the condition This is consistent with our expectations in Section || that G3 has no 

hydrodynamic singularities. 

4.2.2 Shear channel 

According to the discussion in Section ^, equations satisfied by the components ht x , h zx , and 
h ux of the perturbed metric form a closed set. In the radial gauge h uA =0, they read 

H'L = -H' tx + ^H zx + ^-H tx , (4.25b) 
u uf uj 

where Ht x = uht x / (ttT R) 2 , H zx = uh zx j \-kTR) 2 . Using these equations, one finds that 
Zi{u) = <\H tx {u) + voH zx (u) satisfies the following second-order ODE 

„ (w 2 - q 2 f)f - wx> 2 f vx> 2 - q 2 f 
Z x H 7TT7 T Z x H — Z x = . (4.26) 

The exponents of Eq. ( [4.26 ) at u = are and 2, and thus asymptotic behavior of Z\ near 
the boundary is 

Z x {u) = A {1) (1 + ...)+ B {1) u 2 + . . . , (4.27) 

where ellipses denote higher powers of u. The relevant part of the boundary gravitational 
action is given by Eq. (6.19) of [15]. Expressed in terms of the gauge-invariant variable Z\, 
the action is 

5" = — — - — tt— lim / 7^r% , o 7} \ \ Z[( u , k)Zi(u, —k) + contact terms . (4.28) 
8 u^o J (27r) z u{to z — q z /(Mj) 



Proceeding as in Section |4.1| , we find after comparing the expression for the correlator with 

( ggg ) - flOop 

Gi(w, q) = y! . 4.29 

Again, the condition Amicj, q) = is equivalent to Dirichlet boundary condition Z\{u=Q>) = 
0. A typical arrangement of quasinormal frequencies is shown in Fig. ||. 

Hydrodynamic approximation 

In the limit ro < 1, q < 1, the perturbative analytical solution to Eq. ( 4.26| ) satisfying the 
incoming wave boundary condition at the horizon is 15 

Z x (u) =C l f(u)- lW / 2 + ^ +0(» 2 ,q 2 ,»q)) , (4.30) 

15 The argument of footnote ^ regarding the scaling of to, q also applies here. The boundary condition 
constrains vo and q in such a way that the initial assumption that t» and q are of the same order is invalid. 
Using the correct scaling we find that the result (f.32) remains unchanged. 
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Figure 3: Quasinormal spectrum of gravitational fluctuations in the shear channel, shown in the 
plane of complex to = uj/2ttT, for spatial momentum q = q/2nT = 1. The quasinormal frequencies 
coincide with poles of Gi(u;, q), as explained in the text. As q decreases, all poles stay at a finite 
distance away from the real axis, except for the one marked with a large dot. This pole is purely 
imaginary and approaches the origin in the limit q — > 0. The presence of this special quasinormal 
frequency is a manifestation of the diffusive relaxation of transverse momentum density fluctuations 
in the dual Af = 4 SYM theory. 



where C\ is a normalization constant. Expanding for small u, we find the connection coeffi- 
cients 

.4(1) = 1 + ^ + 0(w 2 , q 2 , »q) , = l{V ° 2 2w + 0(tt> 2 , q 2 , tt>q) . (4.31) 

The Dirichlet condition Zi(u=0) = gives the hydrodynamic quasinormal frequency ro = 
— iq 2 /2 + 0(q 3 ). It is interpreted as the dispersion relation for the shear mode, 

c = -i lr] q 2 + 0(q 3 ) , (4.32) 

where 7^ = 1/AirT. For the function Gi(u, q) in this approximation we find 

_ , . vriV 2 T 3 (u; 2 - q 2 ) 

Gi (w, q ) = \ /A 1/ , 4.33 

4(zw — q z /4-Ki ) 

in agreement with the result obtained earlier in [15]. The quasinormal spectrum for frequen- 
cies beyond the hydrodynamic limit was obtained in [11] using a slightly different approach. 

4.2.3 Sound channel 

According to the discussion in Section |||, equations obeyed by the components of the metric 
H tt = uh tt /f{7rTR) 2 , H tz = uh tz / \ttT R) 2 , H zz = uh zz /{7rTR) 2 , H aa = u{h xx + h yy )/{^TR) 2 
form a closed system of equations (in the radial gauge h uA = 0). These equations are lengthy, 
and we present them in Appendix [A]. Using the equations of motion (AT) - ( |A.4Q one can 
show that the gauge-invariant combination 

Z 2 (u) = AwqH tz + 2n> 2 # 2 , + H aa [q 2 (2 - /) - t» 2 ] + 2q 2 /^ (4.34) 
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obeys the following second-order differential equation: 

7 „ _ 3n, 2 (l + ^ 2 ) + q 2 (2n 2 -3?x 4 -3) 
2 u/(3» 2 +q 2 (u 2 -3)) 2 

3t» 4 + q 4 (3 - 4u 2 + u 4 ) + q 2 (4u 5 - Au 3 + 4u 2 to 2 - 6t» 2 ) „ 

H 2 i — 2T^ — ^ Z 2 = 0. (4.35) 

uf 2 (3vo z + q^(u^ — 3)) 



In the limit u — > this equation coincides with Eq. (4.18) obeyed by a minimally coupled 
massless scalar, and thus the behavior of the solution Z2(u) near the boundary is given by 



Z 2 (u) = 4(2) (1 + •••)+ B(2)U 2 + • • • , (4.36) 

where ellipses denote higher powers of u. Using the equations of motion (|A.1[ ) - (|A.4[ ), the 
relevant part of the on-shell boundary gravitational action (quadratic in fluctuations) can be 
written as 16 

Sf = lim J jgf A(v>, q, u) Z' 2 (u, k)Z 2 (u, -k) + S<§, , (4.37) 

where 

3N?Tr 2 T 4 f(u) 

^^ = a<w-A)" <438) 

(2) 

and the "contact term" part does not contain derivatives of the fluctuations (its boundary 
value is given in Appendix |A|). In order to compute the stress-tensor correlation functions, 
we need to solve the "wave equation" (4.35) for Z 2 {u) subject to the boundary condition 

Z 2 ( U =0) = 4wqfl& + 2w 2 H° zz + H° aa (q 2 - w 2 ) + 2q 2 #° , (4.39) 

substitute the result into the action ( |4.37|) and take the appropriate functional derivatives 17 
with respect to the boundary values of the fields H^ t , H^ z , H zz , H® a - For example, for the 
correlator Gu,tt we have 

6 2 s (2) 

Using the expansion ( 4.36 ), we find that the correlators are given by Eqs. (E"3lD - ( HI) 
with 18 

N 2 ir 2 T A B {2) {uJ,q) , x 

G 2 (w, 9) = — — . { ) \ + contact terms . (4.41) 

4(2) {w,q) 

The problem of computing correlation functions in the dual theory is thus reduced to finding 
the connection coefficients 4(2) and 23(2) of the second order ODE Q4.35|) . Zeroes of 4(2) (>»,<?) 

16 The action involving the relevant components of the metric is written in Appendix In order to find 
A(tD,q,it), we form the difference Sg — ^ dujdq/{2%) 2 A{yo,<\,u) Z' 2 {u,k)Z2{u, —k), then use the equations of 
motion (A.l) - (A. 4) to eliminate all derivatives in the difference except H' tt , and find A(XD,%u) by requiring 
that the coefficient in front of H' tt should vanish. 

17 For correct normalization of the coupling between the boundary gravitational fluctuations and the stress- 
energy tensor, see Eq. (3.18) of [16] . 

18 Terms analytic in t» and q are ignored, even if they are divergent as u — > 0. 
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Figure 4: Quasinormal spectrum of gravitational fluctuations in the sound channel, shown in the 
plane of complex to = uj/2ttT, for spatial momentum q = q/2nT = 1. The quasinormal frequencies 
coincide with the poles of G2(uJ,q), as explained in the text. As q decreases, all poles stay at a finite 
distance away from the real axis, except for the ones marked with large dots, which approach the 
origin as q — > (see Appendix ^ for the corresponding dispersion curves). Such behavior of the lowest 
quasinormal frequencies is a manifestation of oscillatory relaxation of longitudinal momentum density 
(as well as energy density) fluctuations in the dual Af = 4 SYM theory. 



appear as poles of the correlators. Finding the poles is therefore equivalent to solving the 
boundary value problem for the gauge invariant variable Z 2 (u) obeying the incoming wave 
boundary condition at the horizon u=l and Dirichlet condition Z 2 (u=0) = at the boundary. 

For q = 0, eqs. ( 4,18| ), ( 4.26 ), ( [4.35 ) all reduce to the equation for a minimally coupled 
massless scalar at zero spatial momentum, and consequently all have the same quasinormal 
spectrum with the asymptotics ( 4.23[ ). This is expected, in accord with rotation invariance 
in the dual field theory, see Eq. ( 2.26| ). For q ^ 0, the spectrum can be found numerically, as 
explained in Appendix a typical arrangement of quasinormal frequencies is shown in Fig. ||. 
However, for small momenta, the lowest quasinormal frequency can be found analytically. 



Hydrodynamic approximation 



In the hydrodynamic regime ro < 1, q <C 1, Eq. ( |4.35| ) for Z^iu) can be solved perturbatively 
in vo and q. Introducing the book-keeping parameter A, rescaling to — > An), q — > Aq, and 
expanding in A <C 1, to first order in A we find 



Z 2 {u) = C 2 f(u) 



-iW/2 



q 2 (l + « 2 



3o> 2 



iwf{u) 



4q 2 



(4.42) 



where C 2 is a normalization constant. Imposing Dirichlet boundary condition Z2 («=()) = 
gives the lowest (|to| <C 1) overtone of the quasinormal spectrum 



vo 



V3 



+ O (q 3 ) • 



(4.43) 
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Figure 5: Real and imaginary parts of the lowest quasinormal (sound wave) frequency as a function of 
spatial momentum. Light curves correspond to the analytic approximation ( |4.43| ) for small q. Dashed 
line is to = q. As q grows, the dispersion curve enters the region where 9Re (t»)/<9q > 1. This, however, 
happens for q ~ 1, where the corresponding singularity of the energy- momentum correlation function 
no longer has the interpretation of a sound wave. 



In the hole-graphically dual finite temperature quantum field theory, Eq. (|4.43| ) appears as 
a pole in the retarded correlator of stress-energy tensors, and is interpreted as dispersion 
relation for the sound wave mode, 

L0{q) =±v s q-iT s q 2 + 0(q 3 ). (4.44) 

The values for the speed of sound v s = l/\/3 and the attenuation constant T s = l/6irT 
coincide with those found previously in [16]. Expanding the solution ( 4.42j ) near u = and 



comparing with Q4.ig| ), we identify the coefficients Am = {<\ 2 — 3» 2 — 2ztoq 2 )/4q 2 , £>( 2 ) = 
(2q 2 + 5itoq 2 — 3ia» 3 )/8q 2 . Thus to leading order in the hydrodynamic approximation we 
obtain 

G2(w,q)= • (4-45) 

The dispersion relation to = to(q) for the sound wave frequency can be determined numerically 
(as explained in Appendix P), and is shown in Fig. [|. 



5. Discussion 



In this paper, we proposed a general approach for identifying quasinormal spectra of asymp- 
totically AdS spacetimes with the poles of the retarded correlators in the holographically dual 
finite temperature field theory. Our demonstration in Section |] that quasi-normal spectrum of 
gauge-invariant perturbations in asymptotically AdS spacetimes has a precise interpretation 
in dual field theory was specific to the case of five-dimensional AdS-Schwarzschild background 
with a plane-symmetric event horizon. Although we expect the same to be true in a more 
general setting as discussed in Section ^ (such as black holes rather than branes, or more 
complicated backgrounds), we did not give a general proof. Explicit computations for other 
backgrounds can be done along the lines of Section ^[ 
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Prom the dual field theory point of view, a noteworthy observation is that the only 
singularities of thermal Green's functions identified in the supergravity approach are simple 
poles. It may indeed be true that real-time thermal correlation functions of gauge-invariant 
operators are meromorphic in complex frequency plane; we leave the analysis of analytic 
structure for future investigation. Of course, this simple nature of singularities is possible 
only at infinite N c ; for example, it is known that low-energy correlation functions of conserved 
currents develop branch cuts in the complex frequency plane, with discontinuities across the 
cuts suppressed in the N c ^cc limit (such cuts are not visible in classical supergravity and 
reflect quantum modifications of the quasinormal spectrum) [28]. In a sense, the situation is 
similar to zero-temperature spectrum of confining gauge theories, when resonances become 
stable as N c ^oo [29]. However, at finite temperature, the poles of the real-time correlation 
functions can not be automatically interpreted as quasi-particles propagating in thermal bath 
- a definite interpretation can be given only after the full spectral density (not just the poles) 
is known. In infinitely strongly coupled N = 4 SYM theory, the quasiparticle interpretation 
is unlikely because of the unique energy scale (temperature). On the other hand, it is possible 
that a' corrections to the quasinormal spectrum will reveal new poles corresponding to heavy 
excitations (whose mass scales with the 't Hooft coupling of the dual field theory) which can 
be interpreted as quasiparticles. 
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A. Equations of motion and boundary action for the sound channel 

Fluctuations of the sound wave mode satisfy a system of differential equations 



Hi 



tz 



2wq (H aa -H U )+ {U ' £ : (» H aa + mH„ + 2q H tz ) 



/ + 2 
3w/ 



q(/ + 2) 
n>(u 3 -3m- 2w 2 ) 
q 2 / 2 (/ + 2) 



<!/(/ + 2) 



(n>-ffaa + + 2q#, 



3«> 2 -q 2 (/ + 2) . 



tz) 



d 2 (/ + 2) 



2«> 2 



/(/ + 2) 
2 



3u - 2i» 2 



/(/ + 2) 



,., f . + . [t» 2 tf aa + ro 2 i^ + q 2 / - H aa ) + 2wqff te ] + ^| 

1 r 9 ^ n 3i» 2 -q 2 (/ + 2) . 
[2w 2 H aa + (u 3 - 3u - 2r» 2 )^] KJ J ' H' H 



+ 



/(/ + 2) 
ro(u 3 — 3u — 2a» 2 
q 2 / 2 (/ + 2) 

1 



q 2 (/ + 2) 



(»F aa + toH zz + 2q# 4z ) , 



(A.l) 



(A.2) 



(A.3) 



^ 2u 2 / 2 (/ + 2) 



6u/(l + n 2 )^ + 2u/ [2t» 2 + q 2 (l + u 2 )] H aa + WufH zz 
8voqufH tz + Auf\ 2 H tt 1 . (A.4) 



The part of the boundary action quadratic in fluctuations is [16] 

1 



B u^o 8 



{3H tt — 12H tz + 2H tt Hn + 2H zz H aa — H zz ) 



2u 



H t z + -7 H l a - H tt Ha + H zz H a 



(A.5) 



where prime denotes the derivative with respect to u, and expressions such as H 2 t are to be 
understood as Htt(u,k)Htt(u,—k). The boundary value of the "contact term" part of the 
gravitational action ( [4.37| ) is given by 

q 2 + 3t» 2 



Sg(0) = lim4 2 )=-^^ 
29q 4 - 30» 2 q 2 + 9» 4 ) 



4(q 2 - to 2 ) 2 



(q 2 — to 2 ) 2 



H%H\ 



4a»q 



/I/O \2 

y^aa) nl^l ^1\ ±± &0' ±± tt .—2 1± aa 1± tz 

2{q z — vo z ) q z — vo z 



tt^ ± tz 



3q 2 + vo 2 rr0 ttQ 9q 4 + 2» 2 q 2 - 3» 4 



rrO ttO 

2\ aa n zz 



+ 



2(q 2 - ro 2 ) 2 

rO \2 



2(q 2 - to 

3q 4 - 18>» 2 q 2 + 7a> 4 2 _ 3q 4 + 14«> 2 q 2 - 9»> 4 , 

\**aa) in 9\9 \^tz 



o 4roq(3q 2 - ro 2 ) 
H tt H zz — — — ti+„ti* 



(q 2 — 1» 

\2 



2\2 1± tz ±± zz 



4(q 2 -w 2 ) 2 



(q 2 - to 2 ) 2 



(A.6) 
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Equation ( A. 6 ) appears to contain more than just contact terms, since there is a pole at 
|t»| = q. This pole, however, is artificial — it reflects the normalization of Z 2 , and cancels in 
the final expression for the correlators. 

(2) 

If one chooses to keep track of contact terms in field theory correlators, one can use Sqj, 
to reproduce contact terms computed earlier in [16]. For example, 

c - 2qA C (u a) 1 ^r(Q) (A7) 
G tm - 3(g2 _ w2)2 G 2 ( W , q) - 4 {6Rg)2 (A.7) 

in the hydrodynamic limit becomes 

3iV c W(3^-5 g 2 ) 
~ 8(^ _ 3^) ' 

which coincides with the result of [16]. 
B. Frobenius solution 

To find the full quasinormal spectrum, one has to analyze wave equations ( |4.5| ), ( (4.18| ), ( [4.26 ), 



19 



and ( 4.35 ) in more detail. All these equations are Fuchsian ODEs with k singular points 
two of which correspond, respectively, to the horizon (u = 1) and the boundary (u = 0). For 
all of the above equations, the exponents at the horizon are equal to ±iw/2, corresponding to 
two local solutions representing outgoing (incoming) waves. The solution obeying incoming 
wave boundary condition at the horizon can be represented as a power series around u=l, 

00 

Z(u) = (1 - u)- iv °l\l + u)-™' 2 «n(»,q)(l - u) n , (B.l) 

n=0 

where Z(u) stands for either of E a , E z , Z±, Z 2 , Z3. The coefficients a n of the series expansion 
obey (k— l)-term recursion relations which can be found by substituting flB.l|) in the original 
differential equations. Quasinormal spectrum is determined by imposing Dirichlet boundary 
condition at u = 0, 

00 

Z(0) =j^a n (»,q) = (B.2) 

n=0 

and solving Eq. ( |B.2j ) numerically taking a sufficiently large but finite number of terms in the 
sum. 20 The spectra of all of the above wave equations are qualitatively similar, except for the 
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For the transverse Maxwell equation (1.5a) satisfied by E a , and for the scalar channel wave equation 



( 4. IS ) satisfied by Z3, the number of singular points is k=4, corresponding to u = 0, ±l,oo. For the lon- 



gitudinal Maxwell equation (4.5b) satisfied by E z , and the shear channel wave equation (4.26) satisfied by 



Zi, the number of singular points is k=6, corresponding to u = 0, ±1, ± v/l — t0 2 /Cf 2 , 00. For the sound 



channel wave equation (4.35) satisfied by Z3, the number of singular points is also k=6, corresponding to 

u ■■ 



20/ 



0, ±1, ±^/3(l -t» 2 /q 2 ),oo. 

Careful readers may note that the expansion (B.l) is guaranteed to converge only inside a circle of radius 



p around u = 1 (in the complex u plane), where p is the distance to the nearest singular point, which may 
become less than one (distance to the boundary) for some values of to and q. However, even for such values of 
to and q, our numerical results for diffusive and shear modes are in agreement with previous calculations [11] 



where this issue did not arise, thus suggesting wider applicability of the expansion (B.l) 
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Figure 6: Real and imaginary parts of three lowest quasinormal frequencies as function of spatial 
momentum. The curves for which ro— >0 as q— >0 correspond to hydrodynamic sound mode in the dual 
finite temperature 7V=4 SYM theory. 



behavior of the lowest (hydrodynamic) frequency which is absent for E a and Z3. For E z and 
Z\ , hydrodynamic frequencies are purely imaginary (given by Eqs. ( gig ) and ( p~32|) for small 
u) and q), and presumably move off to infinity as q becomes large. For Z2, the hydrodynamic 
frequency has both real and imaginary parts (given by Eq. ( [4.44 ) for small u and q), and 
eventually (for large q) becomes indistinguishable in the tower of other eigenfrequencies. As an 
example, dispersion relations for the three lowest quasinormal frequencies in the sound channel 
(including the one of the sound wave) are shown in Fig. ^. The tables below give numerical 
values of quasinormal frequencies for q = 1. Only non- hydrodynamic frequencies are shown 
in the tables. The position of hydrodynamic frequencies at q = 1 is ro = — 3.250637i for the 
.R-charge diffusive mode, ro = — 0.598066i for the shear mode, and ro = ±0.741420 — 0.286280i 
for the sound mode. The numerical values of the lowest five (non-hydrodynamic) quasinormal 
frequencies for electromagnetic perturbations are: 









Transverse channel 


Diffusive channel 








n 


Rero 


Imro 


Rero 


Imro 








1 


±1.547187 


-0.849723 


±1.147831 


-0.559204 








2 


±2.398903 


-1.874343 


±1.910006 


-1.758065 








3 


±3.323229 


-2.894901 


±2.903293 


-2.891681 








4 


±4.276431 


-3.909583 


±3.928555 


-3.943386 








5 


±5.244062 


-4.920336 


±4.946818 


-4.965186 




or g 


ravitational perturbations 


are: 










Scalar channel 


Shear channel 


Sound channel 


n 


Rero 


Imro 


Rero 


Imro 


Rero 


Imro 


1 


±1.954331 


-1.267327 


±1.759116 


-1.291594 


±1.733511 


-1.343008 


2 


±2.880263 


-2.297957 


±2.733081 


-2.330405 


±2.705540 


-2.357062 


3 


±3.836632 


-3.314907 


±3.715933 


-3.345343 


±3.689392 


-3.363863 


4 


±4.807392 


-4.325871 


±4.703643 


-4.353487 


±4.678736 


-4.367981 


5 


±5.786182 


-5.333622 


±5.694472 


-5.358205 


±5.671091 


-5.370784 
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